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Long standing problem in out-of-equilibrium thermal 
field theories are pinching singularities. We find that the 
expressions suspect of pinching require loop particles to be 
on the mass shell. This fact, with the help of threshold ef- 
fect and similar effect due to spin, leads to the elimination 
of pinching in single self-energy insertion approximation 
to propagator in all propagators appearing in QED and 
QCD under very mild restrictions on particle densities. 

This, together with the cancellation of collinear singu- 
larities, allows the extraction of useful physical informa- 
tion contained in the imaginary parts of the two loop dia- 
grams. 

In some cases of interest (it — p interaction, electro- weak 
interaction, decay of Higgs particle, ...) none of the men- 
tioned mechanisms works and one has to resort to the 
resummed Schwinger-Dyson series. These cases are more 
sensitive to the limitations related to the finite time range. 



I. INTRODUCTION 

Out of equilibrium thermal field theories have re- 
cently attracted much interest. From the experimen- 
tal point of view, various aspects of heavy-ion colli- 
sions and the related hot QCD plasma are of con- 
siderable interest, in particular the supposedly gluon- 
dominated stage. 



'Talk given at the 5th International Workshop "Ther- 
mal Field Theories and their Applications", August 10-14, 
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Contrary to the equilibrium case^ I!JEf where pinch, 
collinear, and infrared problems have been -success- 
fully controlled, out of equilibrium theory' Wa^ has 
suffered from them to these days. However, progress 
has been made in this field, too. 

Weldon has observed that the out of equilibrium 
pinch singularity does not cancel; hence it spoils an- 
alyticity and causality. The problem gets worse with 
more than one self-energy insertions. 

Bedaque has argued that in out of equilibrium the- 
ory the time extension should be finite. Thus, the 
time integration limits from -co to +oo, which are 
responsible for the appearance of pinches, have to be 
abandoned as unphysical 0. Similar argument, refer- 
ring to the Fexmi's "golden rule" is given by Greiner 
and Leupold B . „ 

Le Bellac and Mabilat have shown that pinching 
singularity gives a contribution of order g 2 Sn, where 
Sn is a deviation from equilibrium. They have also 
found that collinear singularities canceLm scalar the- 
ory, and in QCD using physical gauges 'liil. but not in 
the case of covariant gauges. Niegawa EJ has found 
that the pinch- like term contains a divergent part that 
cancels collinear singularities in the covariant gauge. 

Altherr and Seibert have found that in massive g 2 (f> 3 
theory pinch singularity— does not occur owing to the 
kinematical constraint V3 . 

Altherr has suggested a regularization method in 
which the propagator is modified by the width 7 which 
is an arbitrary function of momentum to be calculated 
in a self-consistent way. In g 2 (j) A theory, for small de- 
viations from equilibrium, 7 waSrfound to be just the 
usual equilibrium damping rate . 

This recipe has been justified in the resummed 
SchwingeipJ3iison series in various problems with 
pinching oii§. _ , 

Baier, Dirks, and Redlich LJ have calculated the 
7r — p self-energy contribution to the pion propaga- 
tor, regulating pinch contributions by tho-|(iamping 
rate. In subsequent papers with Schiff they 
have calculated .tho-quark propagator within the HTL 
approximation E2rE3 ; in the resummed Schwinger- 
Dyson series, the pinch is naturally regulated by 
ImE R . r-, 

Carrington, Defu, and Thoma til have found that 
no pinch singularities appear in the HTL approxima- 
tion to the cesummed photon propagator . 

Niegawa H has introduced the notion of renormal- 
ized particle-number density. He has found that, in 
the appropriately redefined calculation scheme, the 
amplitudes and reaction rates are free from pinch sin- 
gularities. 

By pinching singularity we understand the contour 
passing between two infinitely close poles: 

/dx 
(x + ie)(x-ie)' 



When e tends to zero, the integration path is 
"pinched" between the two poles, and the expression 
is ill-defined. Integration gives an e _1 contribution 
plus regular terms. Decomposition of (x ± ie)^ 1 into 
PP(l/x) =F iirS(x), gives the related ill-defined 8 2 ex- 
pression. 



The following expression, which is similar to (LI) 
corresponds to the resummed Schwinger-Dyson series 



/ 



dx 



oj(x) 



(x — Y,r(x) + ie)(x — S|j(x) — ie) ' 



(1.2) 



where lu(x) and ui[x) (which appears in (L3)) are, re- 
spectively, proportional to £l(x) and Q(x) where Q(x), 
T.ji(x), and Cl(x) are the components of the self-energy 
matrix. _ _ 

In expression (|L2|), pinching is absent Mia if 
/mSR(i„) =/= at a value of x satisfying x — 
ReTj R (x ) = 0. 

The expression corresponding to the single self- 
energy insertion approximation to the propagator is 
similar to (L2): 



dx- 



uj{x) 



(x + ie)(x — ie) 
One can rewrite the integral as 



dx 

T 

If it happens that 



1 



1 



ui(x) 



hm = K < oo, 

x— >0 X 



(1.3) 



(1.4) 



(1.5) 



then the integral ( |l.4|) decomposes into two pieces 
that, although possibly divergent, do not suffer from 
pinching. 

There are two cases in which the function u)(x) is 
even identically zero in the vicinity of the x = point: 
in thermal equilibrium, because of detailed balance re- 
lations; in massive g 2 (j> 3 theory ojttt of equilibrium, ow- 
ing to the mass shell condition O. The latter mech- 
anism also works in out of equilibrium QED if a small 
photon mass to 7 is introduced. However, this elimi- 
nation of pinching can be misleading: the domain of 
x, where uj(x) = 0, shrinks to a point as to 7 — > 0. We 
shall show that the elimination of pinching also occurs 



in the to. 



cas. 



In this paper Ej we identify two mechanisms lead- 
ing to relation (1.5). They are based on the observa- 



tion that in the pinch-like contribution loop particles 
have to be on mass shell. 

The first mechanism is effective in out of equi- 
librium QED: in the pinch- like contribution to the 
electron propagator, phase space vanishes linearly as 
x — * . In the pinch-like contribution to the photon 



propagator, the domain of integration is shifted to in- 
finity as x — > 0. For distributions disappearing rapidly 
enough at large energies, the contribution again van- 
ishes linearly in the x — > limit. This mechanism is 
also valid in QCD in the cases with massive quarks. 

In out of equilibrium massless QCD, phase space 
does not vanish, but there is an alternative mecha- 
nism: th e sp inor /tensor structure in all cases leads to 
relation (1.5). 

In a few cases, none of the mentioned mechanisms 
works and one has to sum the Schwinger-Dyson series. 
This is the case of the it — p loop in the it self-energy . 
Even in the limit of zer o pio n mass, uj(x) vanishes only 
as \x\ x / 2 and relation (1.5) is not fulfilled. A similar 
problem appears in electroweak interactions involving 
decays of Z and W bosons, decay of Higgs particles, 
etc. Another important case is massless g 2 c/> 3 theory. 
In contrast to massless QCD, massless <? 2 </> 3 theory 
contains no spin factor to provide a q 2 factor necessary 
to obtain (|J). 

The densities are restricted only mildly: they should 
be cut off at high energies, at least as |/co| -3 ' 5 , hi or- 
der to obtain a finite total particle density; for nonzero 
k a , they should be finite; for k near zero, they should 
not diverge more rapidly than |fc | _1 , the electron 
(positron) distribution should have a finite derivative. 

Furthermore, we were unable to eliminate pinches 
related to the double, triple, etc., self-energy insertion 
contributions to the propagator. 

The resummed Schwinger-Dyson series must be free 
from pinching as " for any system which moves towards 
thermal equilibrium and thus behaves dissipatively, 
the full propagator must have some finite width" O . 



II. PROPAGATORS AND THE 
SCHWINGER-DYSON EQUATION 

We start B>0 by defining out of equilibrium ther- 
mal propagators for bosons, in the case when we 
can ignorCj-jijhc variations of slow variables in Wigner 
functions 



k 2 — m 2 + 2ie\k a 



+ 2ir smh 2 96{k 2 - m 2 ), (III) 



D 12 {k) = 27r<5(fc 2 - to 2 ) 

(cosh 2 60(k o ) + sinh 2 0Q(-k o )), {11.2) 
D 2 i(k) = 2ir8{k 2 — to 2 ) 

(cosh 2 9e(-k ) + sinh 2 9Q(k )). (11.3) 



2 



The propagator satisfies the important condition 

= D n -D u -D 21 +D 22 - (II.4) 

To obtain the corresponding relations for fermions, we 
only need to make the substitution 



sinh 2 6»(fc ) -» -sin 2 6»(fc c ). 

In the case of equilibrium, we have 

2 1 
sinh 0(k o ) = n R (ko) 



exp (3\k Q \ — I ' 



(II.5) 



(II.6) 



and similarly for fermions. Out of equilibrium, ns(k ) 
and np(k ) will be some given functions of k . 
We transform to the Keldysh components 



D R (k) = -D 11 +D 21 = 



k 2 — to 2 + 2iek a '' 



D A {k) = —D* R (k) = D R (-k), 



(11.7) 



(Hi 



D K (k) = Du + D 22 = h(k )(D R - D A ) 

= 2n5(k 2 - m 2 )(l + 2 sinh 2 9), 

h(k Q ) = -e(fc )(l + 2 sinh 2 9). (11.9) 
Again for fermions 

D K (k) = 2nS(k 2 - m 2 )(l - 2 sin 2 6). (11.10) 
The proper self-energy satisfies the condition 

= En +Si 2 + E 2i + E 22 . (11.11) 
It is also transformed into the Keldysh form: 

£ R = -(En + E 2 i), E A = S fl , (11.12) 



n = En + s 



22- 



(11.13) 

The "cutting rules" (refs. see also ref. §§] for 

application of the rules out of equilibrium) will con- 
vince us that only on-shell loop-particle momenta con- 
tribute to ImYi R and f2. 

The Schwinger-Dyson equation 

g = G + iGzg, (ii.i4) 

can be written in terms of Keldysh components as 

Qr = G R + iG R V R g R , (11.15) 



Qk = G 



-i (Ga^Gr + Gk^rGr + Ga^aGk) 



(11.16) 



By expanding (11.15), we deduce the contribution 
from the single self-energy insertion to be of the form 



Gr ~ G R + iG R Y. R G R , 



(11.17) 



which is evidently well defined, and the Keldysh com- 
ponent suspected for pinching: 



Gk «G 



K 



-iGa^G r + iGk^rGr + iGa^aGk- 



(11.18) 



The equation for Gr is simple and the solution is 
straightforward: 



Gr 



1 



Gr 1 - i^R 



(11.19) 



To calculate Gk, we can use the solution ( [1.19 ): 
Gk = Ga (^(<Zo)(G a 1 - G B r ) + ifl) Gr- (11.20) 



The first term in ( |ll.20| ) is not always zero, but it does 
not cont ain pinching singularities! The second term 
in ( 11.20 ) is potentially ill-defined (or pinch-like). The 
pinch-like contribution appears only in this equation; 
thus it is the key to the whole problem of pinch sin- 
gularities. In the one-loop approximation, it requires 
loop particles to be on mass shell. 



We start with (|L18|). After substituting ( gL9j ) into 
(11.18), we obtain the regular term plus the pinch-like 
contribution: 



Gk ~ Gkt + Gk p , 



(11.21) 



Gkt = Klo) 

(G R -G a + iGr^rGr - iG A Z A G A ) , (H.22) 

G Kp = iG A nG R , n = n-% )(Efl-E A ). (11.23) 

For equilibrium dens ities, we have E 2 i = e^^S^ , 
and expression ( 11.23 ) vanishes identically. 

Expression ( II.23| ) is the only one suspected of pinch 
singularities at the sing le self-energy insertion level. 
The function f2 in ( 11.23 ) belongs to the type of func- 
tions characterized by the fact that both loop particles 
have to be on mass shell. It is analyzed in detail in 
Sees. Ill and IV (for threshold effect) and in Sec. V 
(for spin effect). With the help of this analysis we 
show that relation (11.23) transforms into 



K 



Gk p — 



.K(q\q Q ) 
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1 



1 



m 2 + 2ieq Q q 2 — m 2 — 2ieq 



(11.24) 



where K(q 2 ,q a ) is &/{q 2 — m 2 ) multiplied by 
spinor/tensor factors included in the definition of 
Gr a- The finiteness of the limit 



lim K(q 2 , q Q ) = K T (q a ) < oo, 

q z — >m J =f 



is important for cancellation of pinches. The index 
=F indicates that the limiting value m 2 is approached 
from either below or above, and these two values are 
generally different. To isolate the potentially diver- 
gent terms, we express the function K{q 2 , q Q ) in terms 
of functions that are symmetric {K\{q 2 ,q Q )) and an- 
tisymmetric {K 2 {q 2 , q )) around the value q 2 = m : 



K(q 2 ,q o ) = K 1 (q 2 ,q ) 

+e(q 2 -m 2 )K 2 (q 2 iqo ). 
These functions are given by 

K ly2 (q 2 ,q ) 

= l(K(q 2 ,q )±K(2m 2 -q 2 ,q )). 



(11.26) 



(11.27) 



Locally (around the value q 2 = m 2 ), these functions 
are related to the limits K±(q ) by 



K h2 (q 2 ,q ) = -(K + ( qo )±K_( qo )). 



(11.28) 



As a consequence, the right-hand side of expression 
(11.24) behaves locally as 



GK P (q 2 ,q ) 



--(Ki(q ) + e(q 2 - m 2 )K 2 (q )) 



1 



1 



m 2 + 2ieq Q q 2 — m 2 — 2ieq 



(11.29) 



and the term proportional to K 2 is capable of produc- 
ing logarithmic singularity. 



III. THRESHOLD FACTOR 

In this section we analyze the phase space of the 
loop integral with both loop particles on mass shell. 
Special care is devoted to the behavior of this integral 
near thresholds. The expressions are written for all 
particles being bosons, and spins are not specified; 
change to fermions is elementary. 



Now, starting from ( 11.12 ) to ( 11.13 ), we calculate 
n and ImTiR-. 



n = 2£Jto£ii = 2^|- / <hi\ n (/.;,.k 1 -<i„)F. (TIL] ) 



where 



d A k 



(11.25) <*M = —4n 2 6(k 2 - m 2 D )S((k - q) 2 - m 2 s ), (III.2) 



and 



Nn{k Q , k a - q ) = -.5e(fc (fc - q )) 

+ (.5 + sinh 2 6 D (k ))(.5 + sinh 2 8 s (k - q )), (III.3) 



Im^R = »— I d»N R (k , k - q )F, (III.4) 



and 



N R (k ,k - q Q ) 

= sinh 2 0D(k o )e{k o - q ) 

+ sinh 2 e s (k - q )e(-k ) 

+9(-fc )6(fc - q ) - G(k )Q(q - k ). (III.5) 

F is the factor dependent on spin and internal degrees 
of freedom. 

It is useful to define N^{k 0l k Q — q a ) as 

N a = Nn - h(q )N R . (III.6) 

After integrating over S's, one obtains 

d 1 \ko\dk. {1 _ 2)) 

Mo] \k\ 

and expressions for f2 and ImT, R take general form 
1 = J d^N(k ,k - q )F(q,k ,\k\,qk), (III.8) 

where \k\ = {k 2 — m^) 1 / 2 , 

qk=\q\\k\z , (III.9) 



<f + k 2 -(q-k) 2 
2\k\\q] 



(111.10) 



0e(O, 2tt) is the angle between vector k? and x axes. 

Let us start with the q 2 > case. Solution of 8(1 - 
z 2 ) gives the integration limits 
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fc oi,2 = 7rj(lo(q 2 +m 2 D - m|) 



2? 



1 



T^I<7l((9 2 -^)(? 2 -e tr )) 1/2 



<?±tr = |tod ± ms| 



(111.11) 



(III. 12) 



Assume now that q tr ^ 0. In this case, at the thresh- 
old, the limits shrink to the value 



qo(gl 



mf) 



2ql 

We define the coefficient ci by 
1 



Cl 



4k1 



d(/)N(k otr ,k otr —q )F. 



(III. 14) 



Now the expression (frll.8 ) can be approximated by 
1 *s Cl {\k\ 2 - \k\i) 



* Cl (0(q 2 - q 2 +tr ) + e(-q 2 + q 2 _ tr )) 

q ((q 2 -q 2 +tr )(q 2 -q 2 _ tr )y/ 2 
Q 2 



(111.15) 



Relation (III. 15) is the key to further discussion of the 
threshold effect. 

We obtain this also for higher dimension (D=6, for 
example) . 

Owing to QIII12D and ( pEII.15| ), the function 
X(q 2 , mjj, m 2 s ) has the following properties important 
for cancellation of pinches. 

It vanishes between the thresholds, i.e., the domain 
(mo—ms) 2 < q 2 < (mc + ms) 2 is forbidden (X = 0). 
If it happens that the bare mass m 2 belongs to this 
domain, the single self-energy insertion will be free 
of pinching. In this case, multiple (double, triple, 
etc.) self-energy insertions, will also be free of pinch- 
ing. Massive A</> 3 theory 113 is a good example of this 
case. 

It is (in principle) different from zero in the allowed 
domain q 2 < (mjj — ms) 2 and (ran + ms) 2 < q 2 ■ 
In this case, one cannot get rid of pinching. This 
situation appears in the tt — p interaction EJ . 

The behavior at the boundaries (i.e., in the allowed 
region near the threshold) depends on the masses mo 
and ms and there are a few possibilities. 

If both masses are nonzero and different (0 ^ mo ^ 
ms ^ 0), then there are two thresholds and X behaves 
as (q 2 —q 2 .^) 1 / 2 in the allowed region near the thresh- 
old q\ tr - For m 2 = q 2 r , the power 1/2 is not large 
enough to suppress pinching. 

If one of the masses is zero (mo ^ 0, 7715 = or 
m,D — 0, ms ^ 0), then (111.15) gives that the thresh- 
olds are identical (i.e., the forbidden domain shrinks 



to zero) and one obtains the {q 2 



il) 1 



behavior near 



m D . This case (for m = m D ) is promising. The elim- 
ination of pinching in the electron propagator, consid- 
ered in Sec. IV, is one of important examples. 

If the masses are equal but different from zero 
(mo = ms 7^ 0), then there are two thresholds 
with different behavior. The function X behaves as 
(q 2 — q^^) 1 / 2 in the allowed region near the thresh- 
er = 4m 2 D ,£ 
pinching in the supposed case m 2 



old ai,^ = 4r7in,and this behavior cannot eliminate 



Am 2 D 



However, at the other threshold, namely at q_ 



(III. 13) 0, the physical region is determined by q 2 < 
and the above discuss ion do es not apply. In fact, 
the integration limits ( III. 11) are valid, but the re- 
gion between k 1 and k Q 2 is now excluded from 
integration. One has to integrate over the domain 
(—00, k a 1) \J(k 2, +00). This leads to the limitation 
in the high-energy behavior of the density functions. 
An important example of such behavior, elimination 
of pinching in the photon propagator (m 7 ), is dis- 
cussed in Sec. IV. 

If both masses vanish (mjj = ms = 0) , the thresh- 
olds coincide, there is no forbidden region and no 
threshold behavior. The behavior depends on the spin 
of the particles involved. For scalars, the leading term 
in the expansion of I does not vanish. Pinching is not 
eliminated. 

The case of vanishing masses (mp = ms = 0) for 
particles with spin exhibits a peculiar behavior. In all 
studied examples (see Sec.V for details), X behaves 
as q 2 as q 2 — > 0, which promises the elimination of 



pinching. 



IV. PINCH SINGULARITIES IN QED 

A. Pinch Singularities in the Electron Propagator 

In this subsection we apply the results of preced- 
ing section to cancel the pinching singularity ap- 
pearing in a single self-energy insertion approxima- 
tion to the electron propagator. To do so, we have 
to substitute mjj — m, ms = 0, sinh 2 &D(k ) —* 
~n e (k ), sinh 2 Os(fc -g ) — > n 7 (k -q ), and h(k Q ) = 
— e(fc D )(l — 2n e {k )), where n e and n 7 are given non- 
equilibrium distributions of electrons and photons in 
relations ( pL5| ) , ( pL6| ) , and ( pOp . The thresh- 

olds are now identical (q± tr = m 2 ), and the integra- 
tion limits satisfy 



\k\2-\k\^ q - 2 {q 2 -m 2 )). 



(IV.l) 



At threshold the limits shrink to the value k tr 
q a , \k\tr = \q\- 



Then, with the help of (III. 14), we define 



5 



#(<Z 2 ,9o) = 



(q 2 — to 2 ) 

d(f>Nn(k otr) k otr - q Q ) 



167T 2 | (q' 2 — TO 2 ) 

W + m)f(£ + m)(|fc| 2 -|fc|i). (IV.2) 
For q 2 ^ 0, we can decompose the vector k as 



fc = — 5— g h — 35—9 + kt 

g z q* 
, Qo -s + to 2 + 9 2 fc . 

— v — T ' 

where, in the heat-bath frame we have 
q = (q o ,0,0, \q\), q = (|g|, 0, 0, g G ), 



lim (n e (q ) + n e (k ) - 2n e (q )n e (k ) . (IV.8) 

The integration limits imply that the limit fc — •> q 
is taken from below for q 2 > to 2 , and from above 
for q 2 < to 2 . The two limits lead to different val- 
ues of N^(q o ,0). This leads to the discontinuity of 
tf(q 2 ,m 2 ,q ) at the point q 2 



Only the first term in ( IV.8| ) can give rise to prob- 
lems. We 
rewrite it as limjfe ^>o ^2fc rt 7 (fc ) d"e(^°+g°) ^ _ re _ 



conditions: 



(IV. 3) lation (IV. 7) should be valid at any q Q we find two 



lim k n^(k ) < 00, 



(IV.9) 



qq = 0, q 2 = -q 2 . 



(IV.4) 



In calculating the term proportional to (1 — a), where 
a is the gauge parameter, we have to use the trick 



((k-q) 2 ±ter 2 
d 



lim 

m-y — >0 



dm 2 



{{k-qf±ie)-m 2 )- 1 



(IV.5) 



Finally, we obtain The "sandwiched" trace factor ]f 
calculated with loop particles on mass shell: 

(Jl + m)f(Ji + m) = 2m{q 2 + to 2 + 2m^) 



+ (q 2 -m 2 ) 



1 2 
q z - TO 



go (g 2 +TO 2 ) fc 



-(1-a) 



(g 2 - to 2 ) 



2<Z 



(IV.6) 



Now we can study the limit 



Jjf(9„) = lim K(q 2 ,m 2 ,q ) 

q 2 -^m l 



(Jj + to) 



5o 



27r|(f|TO ; 



-Nn(k tr, ko tr - g )- (IV.7) 



It is easy to find that If{q ) is finite provided that 
to 2 ^ and N(i(q o ,0) < 00. The last condition is 
easy to investigate using the limiting procedure: 

N u(Qo,0) = lim N^(k ,k - q Q ) 

= lim 2n n (k - q )(n e (q ) - n e (k )) 

k —*q 

+ lim (n e (q ) - n e (k ) - e(q )e(k - q Q ) 

k Q — >q Q 



, dn e (q ) 
dq a 



< 00. 



(IV.10) 



Under the very reasonable conditions ( IV.E ) and 
(1V.10) the electron propagator is free from pinching. 

It is worth observing that If(q ) is gauge indepen- 
dent, at least within the class of covariant gauges. 



B. Pinch Singularities in the Photon Propagator 

To consider the pinching singularity appearing in 
a single self-energy insertion approximation to the 
photon propagator, we have to make the substitu- 
tions mo = to = ms, sinh 2 0d(£; o ) — » —n e (k ), 
sinh 2 6s(fc - q ) — * -n e (k - q a ), and h(k a ) = 
— e(k )(l + 2n~ l (k )). There are two thresholds, but 
only q\ tr = and the domain where q 2 < are rele- 
vant to a massless photon. The integration limits are 
given by the same expression ( III. 1 1 ) , but now we have 



to integrate over the domain (—00, k a 1) (J(fc 2, +00). 
As q 2 — > —0, we find (k a 1 — > —00) and (k a 2 — > +00). 
The integration domain is still infinite but is shifted 
toward ±00 where one expects that the particle dis- 
tribution vanishes: 



K^(q 2 ,q )^ U w -(l- a ) 



q 2 - 2iq Q e 



-jT [9au-(l-a)- 



q 2 + 2iq a e 



ld>T{ 2 \q\q 2 



k a i />oc 



N^(k , k a - q a ) g w - (1 - a) 



9m 1p 



q 2 - 2iq a e 



F pa 3av - (1 - a) 



q 2 + 2iq a € 



(IV. 11) 



G 



To calculate F^ v for the e — e loop, we parameterize 
the loop momentum k by introducing an intermediary 
variable I perpendicular to q. m is the mass of loop 
particles: 



k = aq + I, q.l = o, k 2 = (k — q) 2 = m 2 



,2 _ 2 22 _k 2 + q 2 -(k-q) 2 



I = m — a q , a — 



2q 2 



(IV. 12) 



After all possible singular denominators are canceled, 
one can set a = 1/2. 

Fg = -Tr{Jt + m)^$-k + m)Y 
' 4m 2 q 2 



-A^(q) 



+ ^ (^(4fc (fc - q ) - 4m 2 - q 2 )A^{q) 
+ ('8(k -^) 2 + 2q 2 )B^(q)^, 



(IV.13) 



For pr ojection operators A, B, C and D see (VII. 3) 
( VII. (j ). Now we obtain 



1 



lbn z \q\q' : 



\k\ 



d(f)N^(k , k -q ) 



-2 ° Apu{q) 



(4k (k - q ) - 4m 2 - q 2 )A^(q) 



q 2 



-(-8(k -^) 2 + 2<f)B^(q) 



(IV. 14) 



In the integration over k a the terms proportional to 
{k 2 q 2 ) n dominate and lim g 2_ +0 \K fiU (q 2 , q a )\ < oo if 



'167T 2 |g1 g 2 \J_ OQ J ko2 J \k\ 
{a + (3k 2 Q q 2 ) J d(j)N a {k , k Q - q a )\ < oo. 

Here N^(k , k — q a ) is given by 
Nn(k , k D - q D ) = -2n e (k D - q Q ) 

(-n y (q ) - n e (k )) - n y (q ) - n e (k ) 

-e(q )e(k - q ){-n 7 (q ) + n e (k ) 

+2n 7 (g )n e (fc )). 



(IV. 15) 



(IV.16) 



Assuming that the distributions obey the inverse- 
power law at large energies n e (k ) oc \k \~ s " and 
ns(k ) oc |fc | -<5 % we find that the terms linear in den- 
sities dominate. Thus, for n = 0, 1, one finds 



lk °\ dko \k \ 2n - s (- q 2 y 



q 2 U-oo + J ko2 ) \k\ 

oc (S - 1 - 2n)- 1 (\q\m) 1+2n - S {-q 2 )^-^' 2 . (IV. 17) 



It follows that (IV. 15) is finite (in fact, it vanishes) if 
5 e ,5 s > 3. Similar analysis for electron propagator at 
q 2 < (thus outside of our analysis of pinch singular- 
ities) leads to <5 7 > 3. This is exactly the condition 



d 3 kn lte ^ s (k ) < oo. 



(IV.18) 



Thus the pinching singularity is canceled in the pho- 
ton propagator under the condition that the electron 
and positron distributions should be such that the to- 
tal number of particles is finite. 

Also, in the photon propagator, the quantity 
lim g 2_, K^ [/ (q 2 ,q ) does not depend on the gauge pa- 
rameter. 

Expression JlV.17 ) is not valid for m = 0. 



V. PINCH SINGULARITIES IN MASSLESS 
QCD 

In this section we consider the case of massless 
QCD. One should observe that the massless quarks 
and gluons are an idealisation eventually appropriate 
at the lowest order. In the nonequilibrium HTL re- 
summation scheme both quarks and gluons acquire 
dynamical mass O . Pinching singularities, related to 
massive quarks, are eliminated by the methods used 
in the preceding section. 

Attention is turned to the spin degrees of freedom, 
i .e., t o the function F of the integrand in (111.1) to 



([II. 8). In the calculation of F it has been anticipated 
that the loop particles have to be on mass shell. In 
this case, F provides an extra q 2 factor which suffices 
for the elimination of pinching singularities. 
The integration limits are now 



k oi,2 = ~(q T \q\) 



(V.l) 



The difference \k\2~ is finite and there is no thresh- 
old effect. 

It is worth observing that for q 2 > 0, we have to 
integrate between k a i and k a 2, whereas for q 2 < 0, the 
integration domain is (—oo,k i)1J(A: 2,+oo). This 
leads to two limits, lim 9 2^ ±0 K {q 2 , q a ) = K±(q ), in 
all cases of massless QCD. 
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By inspection of the final results (|V.3| ),( |V.4| ), and 
( |V.5| ), we find that the case q 2 < requires integra- 
bility of the function k 2 N^(k , k a — q a ) leading to the 
condition (IV.1£) on the quark, gluon, and ghost dis- 
tribution functions. 

The function K^(q 2 , q Q ) related to the gluon prop- 
agator is the sum of the contributions from various 
loops, where the terms in the sum are defined as 

K^(q 2 ,q ) = (g^p - (1 - a)D Rw ) 



q- 



-{9av - (1 - a)D Aau ). 



(V.2) 



The tensor F related to the massless quark- antiquark 
contribution to the gluon self-energy is 



1 <M 



^f{jA - 7..)-7-)-l'"'(</' 



-(-8(* " f ) 2 + 2f)B^(q)\ + 0^(k T )J . (V.3) 

As contains only A and B projectors, the result 
does not de pend on the gauge parameter. 

Relation (VJi) contains only terms proportional to 
q 2 , and \im q 2^ K fiu (q 2 , q Q ) is finite. 

For the ghost-ghost contribution to the gluon self- 
energy, the tensor F is given by 



F 



7 9 n = -^N c k^k^qY 



A^(q) 



. T q- 4k (k -q )+q 2 
= -o ab N c — 

T 

-{ ko -^fB^{q)-^D^UO f(n/,-,i). (V.4i 

The tensor F for the gluon-gluon contribution to 
the gluon self-energy is 

pnv _ 5abN c 
99 2 



(g^iq + k) T - g ar (2k - q)» + g Tf *(k - 2q) IT ) 
(k - q) a {k - q) p 



gap - (1 - a) 



(k - q) 2 ± 2i(k - q D )e 



(g vp (q + ky - g m (2k - qf + g^(k - 2q)") 

g Tr ,-(l-a) 
S ab N c q 2 



k 2 ±2ik a e / 

, 7:M (i0(fc o -|) 2 + ^ 2 )^(<z) 



-(-10(fc o - |) 2 + if)B^(q) - ^D^(q) 



-(l-a)[ -A"" - - ^C" 

• 2 \q\ 



2 2 

(1 - a) 2 ( -%A» V + 2%B» V - 2%C> 1V 
• q- q- |g| 



2D^ + 0^{k T ] 



(V.5) 



Expressions (V.3), (V.4), and (V.5) for the ghost- 
ghost, quark- antiquark, and gluon-gluon contribu- 
tions to the gluon self-energy contain only terms pro- 
portional to q 2 . The function K^(q 2 ,q ) approaches 
the finite value K /Jlt/ (±, q ). 

Thus we have shown that the single self-energy con- 
tribution to the gluon propa gator is free from pinching 
under the condition ( IV. 18 ) . 

The K spinor for the quark-gluon contribution to 
the massless quark propagator is defined as 



If(q 2 ,qo)=A 



(V.6) 



In the self-energy of a massless quark coupled to a 
gluon the " sandwiched" spin factor frjffa is given by (as 
the term proportional to fc T vanishes after integration, 
we drop it). 



g^v 



N 2 



1 



2N C 

(1 - a)(k - q)p(k - q) v 



(k - q) 2 ± 2i{k Q - q a )e 



N~ 



1 



2N n 



-fy-^+2^- 1 —^(-fy+^) 

m m 2 \q\ 



(V.7) 



which contains the d ampi ng fa ctor q 2 . 

By inserting ( |V.7|) into ( V.6 ), we obtain ( II.25| ) free 
from pinches. 

To calculate rf{q ), we need the limit 



lim 



q 2 ^0 q z 



. N 2 - 1 2{k - q ) 

°ab—ZT- T fl- 



2N C 



<h 



(V.8) 



From (|V.q ) we conclude that rf{q ) does not depend 
on the gauge parameter. 

Omitting details, we observe that pinching is ab- 
sent from the quark propag ator, also in the Coulomb 
gauge, with the same limit (V.8). 

The K factor for the ghost-gluon self-energy contri- 
bution to the ghost propagator is defined as 
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K(q 2 ,q ) = -o- 



(V.9) 



The F factor for the ghost-gluon contribution is 

Fg h g = 5 ab N c Wq" 

(1 - q)(fc M - qp){k v - q v )" 



9 /J.* 



S ab N c 



(k - q) 2 ± 2i(k - q )e 
„2 



(V.10) 



The factor q 2 ensures the absence of pinch singularity 
and a well-defined perturbative result. 

The K factor for the scalar-photon self-energy con- 
tribution to the scalar propagator is defined as 



K{q 2 ,q ) 



Q 



(V.ll) 



The F factor for the massless scalar-photon contribu- 
tion to the scalar self-energy, 

F S7 = (q + kr(q + ky 

clearly exhibits the q 2 damping factor! 

VI. CONCLUSION 



Studying the out of equilibrium Schwinger-Dyson 
equation, we have found that ill-defined pinch-like ex- 
pressions appear exclusively in the Keldysh compo- 
nent (Gk) of the resummed propagator (11.20), or in 
the single self-energy insertion approximation to it 
( 1I.23| ). This component does not vanish only in the 
expressions with the Keldysh component ( 11.13 ) (f2 or 
f2 for the single self-energy approximation) of the self- 
energy matrix. This then requires that loop particles 
be on mass shell. This is the crucial point to eliminate 
pinch singularities. 

We have identified two basic mechanisms for the 
elimination of pinching: the threshold and the spin 
effects. 

For a massive electron and a massless photon (or 
quark and gluon) it is the threshold effect in the phase 
space integration that produces, respectively, the crit- 



ical q 2 



or q damping factors. 



In the case of a massless quark, ghost, and gluon, 
this mechanism fails, but the spinor/tensor structure 
of the self-energy provides an extra q 2 damping factor. 

We have found that, in QED, the pinching singu- 
larities appearing in the single self-energy insertion 
approximation to the electron and the photon prop- 
agators are absent under very reasonable conditions: 



the distribution function should be finite, exception- 
ally the photon distribution is allowed to diverge as 
k^ 1 as k Q — * 0; the derivative of the electron distri- 
bution should be finite; the total density of electrons 
should be finite. 

For QCD, identical conditions are imposed on the 
distribution of massive quarks and the distribution 
of gluons; the distributions of massless quarks and 
ghosts (observe here that in the covariant gauge, the 
ghost distribution is not required to be identically 
zero) should be integrable functions; they are limited 
by the finiteness of the total density. 

In the preceding sections we have shown that all 
pinch-like expressions appearing in QED and QCD 
(with massless and massive quarks!) at the single self- 
energy insertion level do transform into well-defined 
expressions. Many other theories behave in such a 
way. However, there are important exceptions: all 
theories in which lowest-order processes are kinemat- 
ically allowed do not acquire well-defined expressions 
at this level. These are electroweak interactions, pro- 
cesses involving Higgs and two light particles, a p me- 
son and two 7r mesons, Z, W, and other heavy par- 
ticles decaying into a pair of light particles, etc. The 
second important exception is massless <? 2 </> 3 theory. 
This theory, in contrast to massless QCD, contains 
no spin factors to provide (1.5). In these cases, one 
has to resort to the resummed Schwinger-Dyson se- 
ries. One can also expect that, in these cases, higher 
order contributions become more important and pro- 
vide natural cutoff which reduces the contribution of 
pinch- like terms. In ultimate case this points out to 
the limitations of the method. 

The main result of the present paper is the can- 
cellation of pinching singularities at the single self- 
energy insertion level in QED- anaUQCD-like theories. 
This, together with the reported ErtJ cancellation of 
collinear singularities, allows the extraction of useful 
physical information contained in the imaginary parts 
of the two-loop diagrams. This is not the case with 
three-loop diagrams, because some of them contain 
double self-energy insertions. In this case, one again 
has to resort to the sophistication of resummed prop- 
agators. 



VII. APPENDIX 

We start by defining a heat-bath four-velocity 
Up, normalized to unity, and define the orthogonal 
projector 



9n» - UpUv. 



(VII.l) 



We further define spacelike vectors in the heat-bath 
frame: 



Av^jAv^ — hi — q . 



(VII.2) 
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There are four independent symmetric tensors (we 
distinguish retarded from advanced tensors by the 
usual modification of the it prescription) A, B, and D 
(which are mutually orthogonal projectors), and C: 



(VII.3) 



1 



%q v 



2^q 2 + 2iq a e ' q 2 — 2iq a e 



(VII. 14) 



Br v, u {q) = U IJ ,U v + 



q^qu 

(q 2 + 2iq e) ' 



(VII. 4) 



Cr m „(<?) = 



(-« 2 )!/2 



U.q 

_ iwh, (gp +q 2 )q^qu 

k 2 cpiq 2 + 2iq„e) 



(VII.5) 



Dr ^(q) 



q^qu 



q 2 + 2iq Q e 



In addition to the known multiplication U proper- 
ties (for convenience we drop g-dependence) 



AA = A, B RtA B R>A = B RA , 

Cr.aCr^a = ~(Br.a + Dr,a), 
Dr,aDr,a — Dr,a, 

AB = BA = AC = CA = 0, 
AD = DA = BD = DB = 0, 

(Br,aCr,a)hv = {Cr^aDr^a)hv 
— {Cr^aBr^a)u^ — {Dr^aCr,A)uh 

= q^v 

q 2 ± 2iq € : 
we need mixed products 

Br^aBam = ~^{Br, + Ba), 



Dr A D a ,r = -(Dr + D a ), 



(Br : aCa,r)imv = (C 'r : aD a,r) p.v 



(VII.7) 



(Vlli 



(VII.9) 



(VII.10) 



(VII. 11) 



Cr, a C A ,r = --{B R + B a + D r + D a ), (VH.12) 



(VII. 13) 
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